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The goal of this study is to find a prescription for defining parton distributions (PDFs) 
which are most appropriate for use in those codes where only LO matrix elements 
(MEs) are used, as in many Monte Carlo generators. We describe a modification of LO 
(-H , PDFs, based on using as at NLO, a specific prescription of the coupling in the QCD 

O ■■ evolution, and violation of the momentum sum rule. We compare results with the truth 

- the prediction using NLO for both MEs and PDFs, and the standard LO prediction, 
finding that the modified PDFs generally produce the best results. 
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1 Introduction 



^ ' The talk [T] considers the use of PDFs in LO Monte Carlo generators. It is well known that 

CN I PDFs extracted at different orders of perturbative QCD have large differences in certain 

f^ ' regions of x. This is because missing higher order corrections, both in the parton evolution 

and in the MEs, play an important role in the extraction of the PDFS by comparison to 
experimental data. Traditionally, LO PDFs are supposed to be the best choice for use with 
f^ I LO MEs as implemented in most Monte Carlo programs. However, another viewpoint has 

^TS^ ' recently been put forward, namely it has been suggested that NLO PDFs may be more 

("^ I appropriate even for MEs at LO ^. The main justification for this idea is the claim that 

NLO corrections to MEs are small, and the total cross-section changes due to the differences 
in PDFs. 

j^ ' In this contribution we propose another approach, which combines the advantages of 

H i both the LO and NLO PDFs. We caU the result the LO* PDFs. Here we present only two 

examples of comparisons using all three approaches, but many more examples arc available 
in our previous article ^. However, in this article we additionally introduce another im- 
provement to the modified LO PDFs, namely a change in the scale for the coupling used 
in the evolution of the partons. This makes the PDFs more consistent with the showering 
in Monte Carlo codes and has been inspired by feedback concerning the original LO* PDF 
approximation [1]. 

2 Parton Distributions at Different Orders 

Let us briefly elucidate why the differences between the PDFs at different perturbative 
orders appear. The difference in the gluon PDF is mainly a consequence of quark evolution, 
rather than gluon evolution. The small- a; gluon is determined by dF2/dhiQ^, which is 
related to the Q^ evolution of the quark distributions. The quark- gluon splitting function 
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Figure 1: Parton distribution functions for u-quarks and giuons - absolute values (top) and 
ratio to NLO (bottom) where the uncertainty band for the NLO PDFs is shown. 

Pqg is finite at small x at LO, but has a small-a; divergence at NLO (and further ln(l/x) 
enhancements at higher orders) . so the small x gluon needs to be much bigger at LO in order 
to fit data. Differences of the same nature appear in quark distributions at LO and NLO. 
Most particularly the quark coefficient functions for structure functions in MS scheme have 
ln(l — a;) enhancements at higher orders, and the high-x quarks arc smaller as the order 
increases. Due to the momentum sum rules this is accompanied by a depletion of the quark 
distribution for x ^ 0.01. This depletion leads to a bad global fit at LO, particularly for 
HERA structure function data, which are very sensitive to quark distributions at moderate 
X. In practice the lack of partons at LO is partially compensated by a LO extraction of much 
larger (then at NLO) as(Af^) ~ 0.130. So, the first obvious modification is to use as at 
NLO in a LO fit to parton distributions. The problems caused due to the depletion of partons 
have led to a suggestion [5] that relaxing the momentum sum rule could make LO partons 
rather more like NLO partons where they are normally too small, while allowing the resulting 
partons still to be bigger than NLO where necessary, i.e for the small-a: gluon and high-x 
quarks. We call the modification the LO* PDFs. The approach does improve the quality of 
the LO global fit. The x^ = 3066/2235 for the standard LO fit, and becomes x^ = 2691/2235 
for the modified fit with the same data set as in [6j and using Q;s(Af|) — 0.120 at NLO. 
The momentum carried by input partons goes up to 113%. We analysed the consequences 
of these distributions by comparing to cross-sections using LO and NLO PDFs in LO Monte 
Carlo generators and to full NLO results - presenting many examples [3]. 

Here we introduce an additional modification to the LO* PDFs. We replace the standard 
QCD scale Q^ by Q^ = z(l — z)Q^ in as dependence on the QCD scale in the splitting 
function Pqq (freezing as{Q^) at 0.5 as z —> 1). This automatically takes into account 
logarithmically enhanced terms [7]. The partons obtained using this further modification 
are called the LO** PDFs. As for the LO* PDFS, the output is based on the MRST LO 
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PDFs 0. This improves the quahty of the global fit shghtly further, i.e. x^ = 2640/2235, 
and the resulting value of a5(7\/|) = 0.115, smaller than for the LO* PDFs. In fig. [1] we 
illustrate all four PDF prescriptions: NLO, LO, LO*, and LO**. LO* and LO** gives more 
quarks and especially gluons for x < 0.1 than LO. They provide fewer large-x gluons than 
NLO, but more quarks at large x (as for LO) and more quarks and gluons than NLO for very 
small X. In general, the LO** modification brings little change compared to LO*, mainly 
a small increase in gluons for x ~ 0.01, but in many cross-sections this is countered by the 
smaller value of the coupling. 

As an example of the applica- 
tion of the LO** PDFs we consider 
the M^-boson production and the 
Higgs boson production, both at the 
LHC. For the simulations we used 
MC@NLO [g, CompHEP [TU] and 
FORTRAN HERWIG [U]. We com- 
pare the use of MEs at LO and the 
PDFs at LO/NLO/LO*/LO** to a 
combination of PDFs and MEs at 
NLO (which we call the truth). The 
total cross-sections for the two examples are shown in table [TJ For W production both LO* 
and LO** are closer to the truth than NLO or LO, which is worst, but LO** gives a slightly 
worse answer than LO*. For Higgs production LO** gives the closest number to the truth 
and NLO is worst. The left of Fig. [2] shows two physically important distributions for the 
first example. The LO PDFs show the worst behaviour in the M^-boson pseudo-rapidity, 
reflecting the PDF depletion for moderate x. Our LO*/LO** modifications do not have 
this drawback and imitate the truth much more accurately. The Pt distribution shows that 
we cannot completely simulate the full NLO result with any set of PDFs, though LO* and 
LO** give the closest normalizations. The right of Fig. [5] displays the same distributions for 
our second example. We see that the shapes of all LO approximations are fine, but LO** 
gives a much better normalization. 
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Table 1: Comparison of cross-sections using different 
prescriptions 



3 Conclusions 

We have suggested an optimal set of partons for Monte Carlo codes, which is essentially 
LO but with modifications to make results more NLO-compatible. We call the modification 
LO* and LO** PDFs. They are based on three effects: the use of the NLO QCD couphng, 
relaxing of the momentum sum rule, and in the latter case a change in the scale used for the 
argument of as for high- a; evolution. The resulting PDFs are large where it is required for 
them to compensate for missing higher order corrections, but they are not correspondingly 
depleted elsewhere. We have compared in detail the different PDF approximations combined 
with LO MEs to the truth, i.e. full NLO, for two processes which probe different types of 
PDF, ranges of x and QCD scales. One can conclude that, in general, the results are very 
positive. The LO** and LO* PDFs provide the best description compared to the truth. In 
[3] we saw that this was generally true, especially for s-channel processes, though LO* gave 
a slight overestimate for t-channcl processes. We have confirmed that the LO** PDFs give 
similar results but are in most cases even a little closer to the truth than the LO* PDFs. 
The improvement compared to the LO or NLO PDFs is particularly the case in terms of the 
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Figure 2: Comparison of Z-boson (left plot) and Higgs (right plot) production. The truth 
and all four approximations are shown. 

normalization, but the shape is usually at least as good, and sometimes much better, than 
when using LO or NLO PDFs. It should be stressed that no modification of the PDFs can 
hope to reproduce successfully all the features of genuine NLO corrections. In particular we 
noticed the recurring feature that the high-pT distributions arc underestimated using the 
LO generators, and this can only be corrected by the inclusion of the emission of a relatively 
hard additional parton which occurs in the NLO matrix element correction. However, we 
propose the use of the LO*, or more correctly, the LO** PDFs if only LO MEs are to be 
used. Both LO* and LO** PDFs are now available in the LHAPDF package [13], their 
names are MRST20071omod and MRSTMCal.LHgrid respectively. 
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